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ABSTRACT 

We present the results of HI mapping with the NRAO^ VLA of one of the most metal-deficient blue 
compact dwarf (BCD) galaxies known, SBS 0335-052, with an oxygen abundance of only 1/40 that of 
the Sun. We study the structure and dynamics of the neutral gas in this chemically young object with 
a spatial resolution of 20'.'5 x 15" (~ 5.4 x 3.9 kpc at an assumed distance of 54.3 Mpc), a sensitivity 
at the 2(7 detection level of ~2.0 K or 7.5 x 10^^ cm~^ and a velocity resolution of 21.2 kms^^. We 
detected a large HI complex associated with this object with an overall size of about 66 by 22 kpc and 
elongated in the East-West direction. There are two prominent, slightly resolved peaks visible in the 
integrated HI map, separated in the East-West direction by 22 kpc (84"). The eastern peak is nearly 
coincident with the position of the optical galaxy SBS 0335-052. The western peak is about a factor of 
1.3 brighter in the HI line and is identified with a faint blue compact dwarf galaxy, SBS 0335-052W, 
with rriB — 19.4, and a metallicity close to the lowest values known for BCDs, about 1/50 that of the 
Sun. The radial velocities of both systems are similar, suggesting that the two BCDs SBS 0335-052 and 
SBS 0335-052W constitute a pair of dwarf galaxies embedded in a common HI envelope. Alternatively, 
the BCDs can be the nuclei of two distinct interacting primordial HI clouds. 

The estimated total dynamical mass, assuming the BCDs form a bound system, is larger than ~ 
6 X 10^ M0. This is to be compared to a total gaseous mass Mgas = 2.1 x 10^ Mq, and a total stellar 
mass Mstar < 10® Mq. Hence, the mass of the SBS 0335-052 system is dominated by dark matter. 
Because of the disturbed HI velocity field and the presence of what might be tidal tails at either end of 
the system, we favor the hypothesis of tidal triggering of the star formation in this system. It can be 
due to either the nearby giant galaxy NGC 1376 or the mutual gravitational interaction of the two HI 
clouds. 

"The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under contract with the 
National Science Foundation. 

Subject headings: galaxies: compact - galaxies: dwarf - galaxies: evolution - galaxies: individual 
(SBS 0335-052) - galaxies: ISM - ISM: HI - galaxies: kinematics and dynamics 



1. INTRODUCTION 

Galaxy formation remains a key issue in observational 
cosmology. While much progress has been made in finding 
large populations of galaxies at high (z > 3) redshifts (e.g. 
Steidel et al. 1996; Dey et al. 1998), truly young galax- 
ies in the process of forming, defined as objects which are 
experiencing their first epoch of star formation, remain 
elusive. The spectra of those distant galaxies generally 
indicate the presence of a substantial amount of heavy 
elements, implying previous star formation and metal en- 
richment. 

We adopt here a different approach for studying truly 
^Deceased 1996, September 22. 



young galaxies. Instead of searching for high redshift ob- 
jects, we look for nearby young dwarf galaxies which are 
possibly undergoing their first burst of star formation. As 
pointed out more than a quarter of a century ago by Sar- 
gent & Searle (1970) and Searle & Sargent (1972), some 
low-mass galaxies in the local Universe may approximate 
galaxies in an early stage of their formation. Later, Thuan 
& Martin (1981) identified blue compact dwarf (BCD) 
galaxies as a class of low-luminosity {Mb > -18) extra- 
galactic systems undergoing a strong burst of star forma- 
tion (SF) while exhibiting extremely low heavy element 
abundances. 
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For many years the best candidate for a truly young 
galaxy has been I Zw 18 = Mkn 116, with a metallicity 
Z of about 1/50 Zq (e.g. Lequeux et al. 1979). Stud- 
ies of this galaxy in the last decade suggest indeed that it 
is chemically unevolved. Its extremely low metallicity de- 
rived from the HII region-like spectrum has been confirmed 
by many subsequent studies (see Izotov & Thuan (1998) 
and references therein). Observations in neutral hydro- 
gen (HI) show that it is a very gas-rich system (Lequeux 
& Viallefond 1980; Viallefond et al. 1987; van Zee et al. 
1998a). Also, its HI gas appears to have an extremely low 
metallicity (Kunth et al. 1994), although this last point is 
considered controversial (Pettini & Lipman 1995; van Zee 
et al. 1998a). 

Until the end of the eighties, I Zw 18 remained the only 
BCD with a heavy element abundance as low as Z0/5O. 
Despite efforts by many investigators to search for BCDs 
with equally low or lower metallicities there remained a 
significant gap in metallicity between I Zw 18 and all other 
known BCDs. This gap was partially filled when Izotov 
et al. (1990), using the Russian 6-m telescope, tentatively 
determined the BCD SBS 0335-052 to have a metallicity 
lower than that of I Zw 18. However, later Melnick et 
al. (1992) concluded that Z = 1/41 Zq in SBS 0335-052, 
making it the second most metal-deficient galaxy known 
in the Universe. This BCD is now considered as one of the 
best candidates for a galaxy in formation. 

A dwarf galaxy experiencing its first burst of star for- 
mation should have an extremely low metallicity in its HII 
regions. Such a galaxy should not show any sign of an 
older stellar population produced in previous SF episodes. 
The gas mass would be expected to significantly exceed 
that in stars. Studies of SBS 0335 052 over the last few 
years show that it appears to fulfill all those conditions 
and that it is very probably a truly young galaxy. We 
briefly summarize its main observational properties : 

1. The oxygen abundance in the brightest HII region in 
the galaxy is fomid to be 12 + log(0/H) = 7.30 ± 0.01 
(with variations from 7.15 to 7.33 on small scales), only 
sUghtly higher than that in I Zw 18 ( Melnick et al. 1992; 
Izotov et al. 1997; Izotov et al. 1999). 

2. Blue underlying, extended low-intensity emission is 
detected in SBS 0335 052 on V , R and /-band images. 
The blue {V — I) and {R — /) color distributions suggest 
that a significant contribution to the emission of the ex- 
tended low intensity envelope is due to ionized gas. It is 
found that the observed equivalent width of H/3 emission 
in the extended envelope is 2 — 3 times weaker than the 
value expected in the case of pure gaseous emission. This 
could be explained by the presence of underlying stellar 
emission from coeval A stars (Izotov et al. 1997; Thuan et 
al. 1997). 

3. Prom the optical colors of the low surface-brightness 
component together with evolutionary synthesis models, 
Papaderos et al. (1998) conclude that the age of the un- 
derlying stellar population is less than ~ 100 Myr, with a 
total mass of ~ 10 Mq , comparable with the total mass of 
young blue stars (Izotov et al. 1997). Near- infrared (NIR) 
colors after correction for the gas contribution are consis- 
tent with a stellar population not older than 4 Myr and 
the possible contribution to the NIR light from an evolved 
stellar population in SBS 0335-052 cannot exceed ~ 15% 
(Vanzi et al. 2000). 



4. No evidence for strong, narrow Lya emission is found 
in a UV spectrum obtained with the Hubble Space Tele- 
scope {HST) (Thuan & Izotov 1997). Instead, a damped 
Lya absorption line is observed which implies an HI col- 
umn density iV(HI) = 7 x 10^^ cm~^, the largest ob- 
served so far for a BCD and comparable to the highest 
column densities observed in Lya clouds in the direction 
of quasars. The resonant 01 A1302 line is also detected. 
Assuming that the O I line is unsaturated and that it orig- 
inates in the neutral gas leads to an oxygen abundance in 
the neutral gas of ~ 1/37000 of the solar value, or about 
1/900 of that in the HII regions around the young super- 
star clusters. However, a more consistent interpretation of 
the UV absorption lines of O, Si and S is that these lines 
originate not in the HI but in the HII gas. If this inter- 
pretation is correct, then the neutral gas in this system 
can be pristine and not polluted with heavy elements at 
all (Thuan & Izotov 1997). 

5. Thuan et al. (1995) and Izotov & Thuan (1999) 
have also argued on the basis of the very small dispersion 
of the C/0 and N/0 ratios in extremely metal-deficient 
BCDs with a metallicity less than 1 /20 solar, that in those 
galaxies, C and N are made exclusively by massive (M > 
9 Mq) stars, because intermediate-mass (3 Mq < M < 9 
Mq) stars have not had time to evolve and release their 
nucleosynthesis products. Izotov & Thuan (1999) suggest 
that BCDs with Z < Zq/20 are young, with ages less than 
~ 100 Myr, the main sequence lifetime of a 9 Mq star be- 
ing ~ 40 Myr. Chemical abundances thus also suggest 
that SBS 0335-052, with Z = Zq/40, might be a young 
galaxy. 

In view of the unusual properties of this candidate young 
galaxy we considered it important to measure its neutral 
gas content and study the HI structure and velocity field. 
Measurements of the integrated HI flux of this galaxy in 
the 21-cm line were carried out with the Nangay radio 
telescope (NRT) in 1993 (Thuan et al. 1999a). From that 
data the total HI mass was estimated which indicates that 
this galaxy is gas rich {Mm ~ 10^ Mq) with less than a 
tenth of its baryonic mass in stars. The HI flux appeared 
to be high enough to allow higher resolution mapping with 
the VLA. In this paper we present the results of VLA ob- 
servations of SBS 0335-052 and discuss the properties of 
its HI density distribution and velocity field. 

The paper is organized as follows. In section 2 we de- 
scribe the observations and data reduction. In section 3 
we derive the HI density distribution and velocity field. 
In section 4 we examine the gas and dynamical mass of 
the system and in section 5 we focus on various aspects 
related to the star formation in SBS 0335-052. The paper 
is summarized in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

The blue compact dwarf galaxy SBS 0335 052 was ob- 
served with the VLA in the 21-cm line of neutral hydrogen 
on December 26, 1994 for an on-source time of 163 min in 
the C configuration, and on May 13, 1995 for an on-source 
time of 80 min in the D configuration. The C configura- 
tion, with baselines between 160 to 15000 wavelengths (or 
0.03 to 3 km) was used to achieve high enough angular 
resolution in HI to resolve features comparable in size to 
the optical size of the galaxy (~15"). The D configura- 
tion, with baselines between 160 to 5000 wavelengths (or 
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Table 1 
Observational parameters 



Parameter 


Value 


Value 


Date 


26 Dec 94 


^ r» TV T r\ r- 

13 May 95 


VLA connguration 


C 




No. of antennae 


26 


26 


Time on source [min] 


163 


80 


Central velocity [km s J 


4043 


4043 


TT ir i_ 

Mall-power beam 


20. ox 15 


61 x45 


Brightness temperature [K] 


O A 
Z.U 


O.zi 


corresponding to 1 mJybeam"^ 






Total bandwidth [MHz] 


3.12 


3.12 


Total no. channels 


31 


31 


Channel width [kms^^] 


21.2 


21.2 


Noise in channel maps [mJy beam"-'^] 


0.5 


1.0 


Total bandwidth [MHz] 


1.56 


1.56 


Total no. channels 


63 


63 


Channel width [kms^^] 


5.3 


5.3 


Noise in channel maps [mJybeam"^] 


1.0 


2.0 



0.03 to 1 km), was used to achieve high surface brightness 
sensitivity for gas distributed at larger scales (the 2a noise 
level corresponds to 7.5 x 10^^ cm~^) and to search for 
companion HI clouds, such as those detected near several 
other BCDs by Taylor et al. (1993, 1995, 1996a), Taylor 
(1997). Table I summarizes the observations. 

The observations exploited the 4IF correlator mode. In 
this mode one IF pair was used to take data over a band- 
width of 1.56 MHz, with 63 channels, giving a channel 
width of 24.4 kHz (5.3 kms~^), while the second IF pair 
was used over a 3.12 MHz bandwidth with 31 channels, 
each having a width of 97.7 kHz (21.2 kms^^). The higher 
velocity resolution provides detailed information on the 
galaxy, whereas the lower velocity resolution data cover a 
larger range in velocity space, at higher S/N ratio. 

The data were edited and calibrated with the NRAO 
AiPS package. The C and D configuration observations 
were calibrated separately and inspected. Once satisfied 
with the result, the data were combined and Fourier trans- 
formed. In the end, data cubes were obtained measuring 
512 X 512 pixels, each pixel measuring 5" on a side. The 
resulting synthesized beam was 20'.'5 x 15'.'0. The noise in 
the final cube, at a velocity resolution of 21.2 km s^^, is 
~ 0.5 mjy/beam in a single channel. 

3. HI DISTRIBUTION AND VELOCITY FIELD 

In Figures I and I we show the individual channel maps 
over the velocity range where significant signal is detected. 
The size of the synthesized beam is indicated by the ellipse 
in the lower left hand corner of the first channel map in 
each of the two figures. Fig. |l| displays the HI emission as- 
sociated with SBS 0335-052 which is traced in the channels 
corresponding to a range in velocity from 4000 to 4085 
kms~^. Despite the low signal-to-noise in the channel 



maps, one can appreciate the rather complex kinematics 
of the gas. There is an underlying rotation, with the east- 
ern edge of the HI complex receding. A prominent feature 
in the data cube is a large, almost face-on spiral galaxy, 
NGC 1376. The channel maps corresponding to the ve- 
locities covered by this object are presented in Fig. |^. Its 
velocity ranges from 4064 to 4234 kms"^. The central ra- 
dial velocity, 4162 km s~^ (Rix & Zaritzky 1995), is about 
120kms~^ higher than that of our target galaxy. 

Adopting a Hubble constant Hq = 75kms~-^ Mpc^-^, 
the redshift distance is 54.3 Mpc (Thuan et al. 1997), so 
that 1" corresponds to 263 pc. If NGC 1376 lies at the 
same distance as SBS 0335-052, their projected distance 
of about 9f5 then corresponds to 150 kpc. 

Fig. I is a grey scale presentation of the integrated inten- 
sity of the HI in the target field. It is obtained by creating 
a zeroth moment map of a blanked data cube, i.e., a data 
cube where all signal below the 2(7 level was deleted and in 
which only those regions are preserved which show emis- 
sion in at least two successive channels. Both SBS 0335- 
052 and NGC 1376 are displayed. Notice the elongated 
structure of the SBS 0335-052 HI cloud with two promi- 
nent concentrations separated by about 84" (or 22 kpc). 
These "peaks" are marked respectively by "E" and "W" 
on the map. They are slightly resolved, as indicated by 
gaussian fitting of the brightness distribution of small re- 
gions centered on these peaks, and are situated approx- 
imately symmetrically relative to the cloud center. The 
western peak is about a factor of 1.3 brighter than its 
eastern counterpart. The HI cloud as a whole is quite 
inhomogeneous . 

Both the E and W HI peaks are embedded in a common, 
elongated envelope with a strong bend at the eastern edge. 
The outermost contour in Fig. 3 corresponds to a HI col- 
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Table 2 

Coordinates of reference stars used 
to tie the optical to the vla 
coordinate frame 



Star a(i950.o) '5(i950.o) 



1 


03^34'°59!71 


-05°13'04'.'0 


2 


03'^35'"02!64 


-05°12'05'.'9 


3 


03''35'"16!93 


-05°13'19'.'7 


4 


03''35'°19!23 


-05°11'57'.'7 


5 


03'^35"21?42 


-05°10'29'/8 


6 


03''35'°26!86 


-05°10'50'.'6 



-05 09 



„ -05 09 



-05 09 



• ' ° ' 3958.3 ■ 
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RIGHT ASCENSION (B1950) 

Fig. 1. — Mosaic of channel maps of the low velocity resolution data showing the HI distribution of SBS 0335-052. Contours represent 1 
(2(t), 2, 4 and 6 mjy beam~^ . The ellipse in the lower left hand corner of the panel on the upper left indicates the beam size. 



umn density of 3 x 10^° cm^^ and a size of about 31x7 
kpc. At the 2(7 level above the noise (represented by the 
lightest grey scale tones in Fig. 3), which corresponds to a 
HI column density of 1.9 x 10^" cm~^, the total extent of 
the HI emission is 66 x 22 kpc. Some additional weak HI 
features may be present in the field around SBS 0335-052, 
including an extended object to the NE. However, these 
detections are tentative at best and should be confirmed 
with future higher S/N observations. The structure of the 
HI feature to the NE of the SBS 0335-052 HI complex at 
(5 = -05° 09' 30" is quite unusual. No cases similar to 
this one were found in studies devoted to a search for HI 
companions near HII galaxies (Taylor et al. 1993, 1995) or 
low-surface-brightness (LSB) dwarf galaxies (Taylor et al. 
1996b). 

We should like to emphasize that an HI envelope of this 
size is rare among BCDs and dwarf irregulars. The typical 



size of HI envelopes around BCDs is usually a few kpc (see 
for example the discussion in Viallefond & Thuan (1983) 
and in van Zee et al. (1998b)). The only known HI enve- 
lope of comparable size is that of HI 1225-1-01 (Chengalur 
et al. 1995). 

Figure ^ represents the integrated profile of the HI cloud 
associated with SBS 0335-052. It was calculated by inte- 
grating all the flux in the region around SBS 0335-052, as 
seen in Figure ^. The dashed lines are used to indicate 
separately the integrated profiles of the eastern and west- 
ern halves of the HI cloud. Whereas the integrated profile 
width measured with the VLA is close, within the errors, 
to that measured with the NRT (Thuan et al. 1999a), 
the VLA integrated flux of 2.46 ± 0.18 Jy kms^^ is about 
a factor of two higher than that measured with the NRT 
(1.28 ± 0.22 Jy kms^^). This difference points to a signif- 
icant loss of integrated flux for the NRT due to its narrow 
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RIGHT ASCENSION (B1950) 

Fig. 2. — Mosaic of channel maps of the low velocity resolution data showing the HI distribution of the bright spiral galaxy NGC 1376. 
Contours represent 1 {2(t), 2, 4, 6 and 8 mjy beam~^ . The ellipse in the lower left hand corner of the panel on the upper left indicates the 
beam size. 



horizontal beam (FWHM = 3.'5) as compared to the larger 
extent of up to ~ 4' of the HI cloud in the EW direction. 
After a first order correction for the NRT beam shape, 
the NRT/VLA flux ratio becomes 0.78 ± 0.27, which is 
consistent within the errors. 

Fig. H is an overlay of the HI contours on top of an op- 
tical S-band image from Papaderos et al. (1998). The 
optical astrometry was done with the six stars marked by 
crosses. Their coordinates were obtained from the Hubble 
Space Telescope Guide Star Catalog (Lasker et al. 1990) 
and are given in Table 2. As shown first by Pustilnik et 
al. (1997) the two HI peaks are clearly identified with two 
faint optical galaxies — the eastern one with the object 
SBS 0335-052 itself, and the western peak with a very 
compact, faint dwarf galaxy with tub = 19.4. A compar- 
ison of the optical and HI peak positions will be given in 
section 4. No 20 cm continuum emission was detected at 
the locations of SBS 0335-052, SBS 0335-052W or NGC 
1376 at the 4 a level of 1 mJy. To avoid as much as pos- 
sible confusion, we will use the same nomenclature for the 
HI components as for the optical counterparts, referring to 
the eastern HI component simply as SBS 0335-052 and to 
the western one as SBS 0335-052W. When referring to the 
entire HI complex we will employ the term SBS 0335-052 
system. 

In Fig. 1^ we show an enlarged image of the integrated 
HI flux distribution in the SBS 0335-052 system, the grey 
scale being a linear representation of the integrated HI 
flux, with superimposed isovelocity lines with a contour 
interval of 5 kms~^. The velocity field of the HI gas is 
quite complex. There is a general gradient across the sys- 
tem from 4085 kms~^ at the eastern edge to 4000 kms^^ 
on the western rim. However, there are clear deviations 
from this general trend near the E and W peaks. This 



is perhaps better illustrated in Fig. [7| which is a position- 
velocity (P-V) diagram made along a line through the two 
peaks and based on the low velocity resolution data. If the 
HI complex is treated as a single system, then there is a ve- 
locity gradient across, with a maximum velocity difference 
of ~ 80 kms^^ between the two sides. If it is considered 
however as two distinct HI clouds, then the velocity gra- 
dients are solid-body in nature, with maximum velocity 
differences of ~ 35-40 kms"""^ within each cloud, 
[hbtp] 

4. RESULTS 

As mentioned in the previous section, we derive an in- 
tegrated HI fiux over the entire region around the cloud of 
2.46 Jykms"^. This corresponds to a total HI mass: Mhi 
= 1.68 X 10^ Mq. Since the blue luminosities of the east- 
ern and western galaxies are respectively 7.3 x lO^L© and 
0.82 xIO^Lq, the Mm/ Lb of the SBS 0335-052 system is 
of order 2.1. 

Accounting for a mass fraction of 0.245 for helium (Izo- 
tov & Thuan 1998), the total gas mass in the system is: 
Mgas = 2.1 xlO^M©. 

With a Salpeter initial mass function, and lower and 
upper mass limits of 0.8 and 120 AIq respectively and us- 
ing the models of Schaerer & Vacca (1998), Papaderos et 
al. (1998) have derived a total stellar mass for the star- 
burst and underlying stellar components of 3.1xlO^M0. 
Lowering the lower mass limit to 0.1 Mq would increase 
the total stellar mass to «3.7xlO^Mo. From mid-infrared 
observations, Thuan et al. (1999b) have found that as 
much three-quarters of the current star formation activ- 
ity in SBS 0335-052 can be hidden by dust. Correcting 
for this effect would bring the total stellar mass to 7.3 x 
IO^Mq and 8.8 x IO'^Mq for lower mass cut-offs of 0.8 
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Table 3 

Some observed and derived characteristics of the SBS 0335-052 system 



Parameter 


0335-052W 


0335-052E 


References 


a(i950.o) (opt)'' 


03'^35"09!53 


03^35'"15!17 


4 


<5(1950.0) (opt) 


-05°12'25"2 


-05°12'27'.'6 


4 


B [mag] 


19.37 


17.00 


1 


Lb [108 LbQ] 


0.82 


7.3 


1 


optical angular size ["]^ 


14x14 


23x20 


1 


optical linear size [kpc] 


3.7x3.7 


6x5.3 


1 


Vopt [kms"^] 


40G9±20 


4060±12 


2,3 


Vm [kms-i] 


4017± 5 


4057± 5 


4 


"(1950.0) (HI) 


03'^35'"09?50 


03'^35'"15?21 


4 


<5(1950.0) (HI) 


-05°12'25'.'5 


-05°12'21'.'5 


4 


Peak A^Hi [10^° cm'^] 


9.98 


7.42 


4 


Peak HI mass surface density [M© pc~^] 


8.0 


6.0 


4 


M"b 


-14.3 


-16.7 


1 


HI angular size ["Y 


128x67 


122x83 


4 


HI linear size [kpc] 


33.7x17.6 


32.1x21.8 


4 


Mhi [Mq/IQS] 


8.9 


8.0 


4 


Mhi/Lb \M^^/LbQ\ 


10.8 


1.1 


4 


Mi.jLu \M. /Lis ■] 


2().<S 


2.2 


4 



''Optical coordinates derived from DSS-2 images. 

''Optical size corresponding to the isophotal level of 27.5 mag arcsec"^ in B (Papaderos et al. 1998). 

■^Maximum observed extent of the associated HI gas. 

'^Lbq is the solar blue luminosity corresponding to Mb0=5.48 

References.— 1. Papaderos et al. (1998). 2. Lipovetsky et al. (1999). 3. Izotov et al. (1997). 4. 
This paper. 
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Fig. 3. — Total HI map of SBS 0335—052 and companion galaxy NGC 1376, based on the low velocity resolution data. Contours represent 
1.9 (5(7), 3.8, 7.6 and 11.4 X 10^" cm"'^. The grey scale is a linear representation of HI surface density ranging from (white) to 16 X 10^" 
cm~^. The ellipse in the lower left hand corner indicates the beam size. 



3 

E 



40 



30 



20 



10 



SBS 0335-052 




3800 



3900 



4000 
Velocity (km/s) 



4100 



4200 



Fig. 4. — Integrated profile of HI fiux (in mjy) for the entire HI cloud (solid line) in the system SBS 0335-052, as derived from the channel 
maps displayed in Fig. 1 for the low velocity resolution data. The total fiux, integrated over this profile is 2.46 Jykms"'^. Dashed profiles 
show integrated profiles for the E and W components separately. 



and 0.1 Mq respectively. 

Since SBS 0335-052W is about one order of magnitude 
less luminous and slightly less chemically evolved, its stel- 
lar mass is expected to be about 10 times lower than that 
of the eastern galaxy. Thus the total stellar mass in both 
dwarf galaxies is at most IO^Mq, or no more than 5% of 
the total gas mass. 

The two peaks in the distribution of the integrated HI 
map are partly resolved with the synthesized beam of 
20'.'5 X 15'.'0 (or 5.4 x 3.9 kpc at the adopted distance of 
54.3 Mpc). The position angle of the beam is PA = -37°. 



We used the AiPS task imfit to estimate the properties of 
both HI peaks. The central brightness of the western peak 
is 25.8Kkms~^, and its FWHM size after deconvolution 
with the beam is (26'.'7 ± 1'.'5) x (15'.'6 ± 1'.'5) at PA = 
53° ± 7°. For the eastern peak we get a central brightness 
of 19.2Kkms^^, and a FWHM size after deconvolution of 
(28'.'4 ± 1'.'7) x (12'.'7 ± 4'.'9) at PA = 53° ± 9°. These 
central brightnesses correspond to a column density of HI: 
iVHi(west) = 10.0 x 10^0 cm-2 = 8.0 Mq pc'^, 
iVHi(east) = 7.4 x lO^o cm^^ = 6.0 Mq pc'^. 
The faintest regions of the outer HI disk (roughly cor- 
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Fig. 6. — Isovelocity lines of HI superimposed on the HI surface brightness map. The grey scale is a linear representation of HI surface 
density ranging from (white) to 11 X lO^" cm ^ . 



responding to a 2f7 noise level over two consecutive chan- 
nels) correspond to column densities of 7.5 x 10^^ cm"^ 
or about 0.6 Mq pc~^. Both HI concentrations are elon- 
gated, oriented at the same position angle, roughly SW 
to NE. Their FWHM angular sizes are very similar and 
correspond to a linear extent of about 7.1 x 3.7 kpc. In 
the SW-NE direction the intrinsic (deconvolved) width of 
both HI peaks is about a factor of two larger than that 
of the VLA beam, but in the perpendicular direction it is 
only 0.7 - 0.8 of the beam width. 



The masses of the eastern and western components, sim- 
ply dividing the entire cloud in two halves, as was done in 
Figure ^, leads to gas masses (corrected for He) of: 

Mffas(west) = 1.11 xIO^Mq, 

Mgas(east) = 0.99 xIO^Mq. 

The peak column density, measured for the E-peak, 
7.4 X 10^° cm~^ is a factor of 10 lower than the value 
derived from Lya absorption based on HST observations 
(Thuan & Izotov 1997). This is also a factor 4 less than 
the peak brightness of HI for I Zw 18 as derived from VLA 
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observations by van Zee et al. (1998a) with a synthesized 
beam ~ 5" (corresponding to a hnear resolution of 0.25 
kpc). Both facts seem to indicate that despite the fact 
that some fraction of the brightest parts of the SBS 0335- 
052 HI complex is resolved by our VLA beam, structure 
at scales significantly smaller than the beam must exist 
and that the structure of the neutral gas in the vicinity of 
the region with high current SF rate is clumpy, explaining 
the difference in the peak HI column density derived from 
fitting the damped Lya line profile (Thuan & Izotov 1997) 
and that derived here. Probably, an inhomogeneous gas 
distribution due to HII regions and SNR shells affect the 
structure of the neutral ISM. Such shell structure is al- 
ready hinted at by the optical image obtained by Melnick 
et al. (1992) and confirmed on HST images by Thuan et 
al. (1997). 

The coordinates of the HI peaks, determined with the 
AiPS task IMFIT, as defined by the centers of gaussian com- 
ponents, are presented in Table ||. Their r.m.s. uncertain- 
ties are (f.l8 x 1'.'3 for the eastern HI peak, and 0!11 x 
O'.'S for the western one. 

The optical position of SBS 0335-052W coincides within 
the errors with the position of the corresponding HI peak. 
For SBS 0335-052, while there is no shift in right ascension 
between the optical and HI peaks within the cited errors, 
there exists a difference in declination A5 = +6'!l. This 
is more than twice the combined error in declination and 
seems significant. To check the reality of this offset, we 
have also compared the optical image with a map made 
only with the C-configuration data. The displacement is 
still present. Taken at face value, the shift between the 



optical and HI peaks for SBS 0335-052 is about 1.6 kpc, 
which corresponds to the offsets reported by, e.g., Vialle- 
fond & Thuan (1983). Note, however that van Zee et al. 
(1998a, b), with higher spatial resolution HI maps, find no 
significant offset between the peaks of the HI and optical 
intensity distributions in the BCD I Zw 18. 

Spectroscopy of the newly discovered galaxy at the po- 
sition of the western peak with the 6~m (Pustilnik et al. 
1997), and with the Multiple Mirror and Keck II telescopes 
(Lipovetsky et al. 1999) reveals an emission-line spec- 
trum with a radial velocity close to that of the HI cloud. 
SBS 0335-052W, with an absolute magnitude Mb about 
— 14.3 is, like SBS 0335-052, an extremely metal-deficient 
object. Its metallicity is 1/50 that of the sun, comparable 
to that of I Zw 18. 

The total dynamical mass is difficult to determine as this 
depends crucially on one's interpretation of the SBS 0335- 
052 system. In this paper we will present two strawman 
hypotheses. The first one is based on the assumption that 
the system is in fact one cloud with two condensations. 
The second hypothesis is that we are witnessing the inter- 
action of two systems and that the observed HI shows a 
tidal bridge between both optical galaxies and tidal tails 
stretching towards the northeast and to the west. 

In the case of one giant cloud, we can estimate the total 
gravitational mass by equating the centrifugal and gravi- 
tational forces at the edges of the disk. We take for the 
maximum rotational velocity 40 kms~^, and a radius of 
the disk equal to 16 kpc (Fig. 7). The resulting estimated 
lower limit for the total dynamical mass is: 

Mdyn{R < 16 kpc) = 5.9 X 10^ Mq. 
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In the case of two interacting clouds we estimate, based 
on the position-velocity diagram (Fig. 7) a full width 
to zero velocity range of 65kms~^ for SBS 0335-052 and 
75kms-i for SBS 0335-052W. The linear sizes of the HI 
clouds are at most 12kpc (this is using the full extent 
rather than the FWHM along the major axis of each ob- 
ject). We then find: 

Md3^„(east) = 1.6 X 10^ Mq, 

Mdj;„(west) = 2.2 X 10^ Mg, 
or a total of Mdyn (east+wcst) = 3.8 x 10'' Mq, 64% of 
that derived on the basis of the more extreme hypothe- 
sis of the objects forming one cloud. It should be noted 
that these estimates are strict lower limits as the inclina- 
tion is completely unknown and was assumed to be 90° 
(or edge- on). 

We have shown that in both hypotheses, we need to 
invoke substantial amounts of dark matter in order to rec- 
oncile the observed baryonic content with the dynamical 
masses derived on the basis of the observed kinematics. 
To obtain yet anothca' (lower) limit to the amount of dark 
matter which might be hiding within the SBS 0335-052 
system, we can take the observed (projected) separation 
and radial velocity difference and calculate the mass which 
this implies assuming the two clouds to be in a bound or- 
bit. This works out to be Mtot = 9 x 10^ Mq (see van 
Moorsel 1987 for details). Therefore, both the individual 
clouds, as well as the SBS 0335-052 system require sub- 
stantial amounts of non-visible matter to make up for the 
inferred masses as based on their dynamics. 

5. DISCUSSION 

5.1. SBS 0335-052 in relation to other candidate young 

galaxies 

It is interesting to compare this system with other can- 
didate young galaxies. There are three galaxies in the lit- 
erature with similar characteristics which have been pro- 
posed as possible young objects: I Zw 18, ESO 400— G43 
and HI 1225-1-01. 

The HI in I Zw 18 has a complex structure with three 
distinct HI components (Viallcfond et al. 1987; van Zee 
et al. 1998a), two of them associated with star-forming 
regions: the NW and SE components in the main body 
of the galaxy and component C to the north (e.g. Du- 
four et al. 1996). But unlike SBS 0335-052, where the 
two centers of star formation are well separated spatially, 
the three star-forming regions in I Zw 18 are connected 
causally, star formation self-propagating from component 
C to the SE component (Izotov & Thuan 1998). 

In the case of ESO 400-G43 (Bergvall & Jorsater 1988), 
the two HI components arc well separated by a projected 
distance of ~ 40 kpc, each component containing an opti- 
cal counterpart. Only one of the optical galaxies has been 
studied in detail. It has been described as a possible young 
galaxy, although its metallicity Z = Zq/8, is probably too 
high for it to be forming its first generation of stars (see 
Izotov & Thuan 1999). 

The HI 1225-hOl system (Salzer et al. 1991; Chengalur 
et al. 1995); at a distance of 20 Mpc, presents several sim- 
ilarities to the HI envelope associated with SBS 0335-052. 
It also has two extended components, separated by a pro- 
jected distance of 98 kpc, and with an HI tidal bridge 
linking them. A dwarf HII galaxy 12 kpc in size is located 



within the 40 kpc NE HI component, approximately at the 
position of the HI peak. The total gas mass of the NE HI 
component of ~ 3 x 10^ Mq, is close to that of the HI 
gas in the SBS 0335-052 system. The dynamical mass is 
a factor of 3 larger than the gas mass. Due to a smaller 
inclination angle {i ~ 30°) and better spatial and velocity 
resolution the authors arc able to distinguish bar-like and 
spiral structure in the NE HI component. 

The very low metallicity, Z ~ 1/20 Zq, in the HII re- 
gion associated with the current SF burst and the very 
low mass fraction in the form of stars suggest that the 
HI 1225+01 system probably is in the process of form- 
ing its first generation of stars. One of the arguments by 
Salzer et al. (1991), who claim that the current SF burst 
in this HII galaxy is possibly not the first one, is that there 
is a significant abundance of nitrogen, which the authors 
assume to be produced by intermediate-mass stars. How- 
ever, as mentioned earlier, the finding of a constant N/0 
ratio with a very small dispersion in the most metal-poor 
BCDs, those with Z < Zq/20 (Thuan et al, 1995; Izotov 
& Thuan 1999) strongly suggests primary production of 
nitrogen by massive stars in galaxies as metal-deficient as 
SBS 0335-052, and thus there is no need for a second burst 
of star formation in the optical counterpart of HI 1225-1-01 
to account for the nitrogen abundance. The second HI 
component has no detectable optical emission, but shows 
a velocity field corresponding to that of an inclined rotat- 
ing disk. Thus, Chengalur et al. (1995) suggest that the 
most probable interpretation for this system is a SF burst 
in a primordial HI cloud (the NE component) due to the 
tidal action of another massive HI cloud. 

The HI properties of the SBS 0335-052 system also 
show some similarities with those of II Zw 40 (Brinks & 
Klein 1988; van Zee et al. 1998b). The metallicity of this 
BCD is Zq/6, considerably higher than those of the 
above-mentioned systems, and suggesting that this galaxy 
is more evolved than SBS 0335-052 and not necessarily its 
analog. Nevertheless it is instructive to briefly compare the 
two objects to discuss their different evolutionary states. 
The II Zw 40 system is most likely a merger, showing two 
HI tidal tails of ^ 6.5 and 15 kpc in size. The extent of the 
HI gas and the HI masses of the tails are reminiscent of 
the situation in SBS 0335-052. However, the optical mor- 
phologies of SBS 0335-052 and II Zw 40 differ. Whereas 
II Zw 40 resembles an advanced merger with optical tidal 
tails (van Zee et al. 1998b), SBS 0335 052 exhibits two 
regular-looking still widely separated dwarf galaxies which 
are perhaps in an earlier stage of what could evolve into a 
merger. Over time, SBS 0335-052 might very well develop 
into a system resembling II Zw 40. 

5.2. Possible star formation triggers 

One of the unsolved problems ever since the discovery of 

the first BCDs (Sargent & Searle 1970) concerns the trig- 
gering mechanism for the star formation. And perhaps 
even more baffling, why has star formation not occurred 
earlier in the lifetime of these objects. 

In all likelihood, tidal triggering is probably responsible 
for the current star formation in the SBS 0335-052 sys- 
tem. We shall again consider two hypotheses: 1) the case 
where SBS 0335-052 is one huge HI cloud containing two 
star-forming centers; and 2) the BCDs are the nuclei of 
two distinct interacting HI clouds. 
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In the case of a single sclf-gravitating HI cloud, the 
tidal triggering is probably due to the companion mas- 
sive galaxy NGC 1376 whose HI map is seen in Figure 
8. The resolution of the map is too low to say much 
about the HI morphology. However we can note that the 
HI gas is, as usual, more extended than the optical disk, 
and that there is a conspicuous lack of HI in the nuclear 
region of the galaxy. A lower limit to the distance be- 
tween the galaxy and the SBS 0335-052 system is their 
observed projected distance of 150 kpc. NGC 1376 is an 
Scd galaxy with nipg = 12.8, Mb = —21.0, an optical di- 
ameter D = 2'.0 and an inclination of 21°. Its HI profile 
is characterized by a full width at 20% of peak intensity 
of 179 kms~^ and a heliocentric radial velocity Vhi = 
4162 kms~^. It is a member of a medium density group 
listed as LGG 103 (Garcia 1993) and composed of at least 
14 galaxies. NGC 1376 is situated on its south-western 
outskirts. The closest known group member is at a pro- 
jected distance of about 0.7 Mpc. Rix & Zaritzky (1995) 
note that NGC 1376 has 3 spiral arms, and Elmcgrcen & 
Elmegreen (1987), in their classification of spiral galaxies, 
assign this galaxy to class 2, which means that its arms are 
"fragmented spiral pieces with no regular pattern" . This 
suggests that the spiral has recently experienced a tidal 
disturbance. 

The VLA HI data does point to a small galaxy, compan- 
ion to NGC 1376, that may be at the origin of that tidal 
disturbance. Fig. 8 shows a small irregular HI blob SW 
of the main body of NGC 1376 at a = 03^34'°32.9" and 6 
= -05°13'51". Its detection is weak, but fairly reliable as 
it can be seen in at least two neighboring channel maps, 
at a heliocentric velocity of 4117±5 kms~^. Its integrated 
flux density is 0.17 Jy km s~^, which corresponds to an HI 
mass of 1.2 x 10^ Mq at the distance of 54.3 Mpc. This 
blob is kinematically detached from the regular velocity 
field of the galaxy and can be identified optically on the 
Deep Sky Survey (DSS) with what appears to be a faint 
elongated dwarf companion, not listed in NED and with 
a diameter of ~ 23" (i.e. a linear size of ~ 6 kpc) at a 
projected distance of 27.5 kpc from NGC 1376. The HI 
profile of the dwarf galaxy is two channels wide so that its 
velocity width is ~ 42 km s^^. This gives a lower limit 
(the inclination being unknown) of ^ 2.5 x 10^ Mq for its 
dynamical mass. 

The total dynamical mass of NGC 1376, as estimated 
from the width of the HI profile (inclination corrected to 
Vrot = 250 kms^^) and the radius of the HI distribution, 
is about 5 x 10^^ Mq. Taking into account the results 
by Zaritsky & White (1994) and Zaritsky et al. (1997) 
on distant dwarf satellites of spiral galaxies and masses 
of their dark matter halos, the total gravitational mass of 
NGC 1376, exerting a tidal force on SBS 0335-052, can 
be up to twice that value, or 10^^ Mq. A detailed under- 
standing of the triggering mechanism of a SF burst due 
to tidal forcing such as due to NGC 1376 is still lack- 
ing. It will require a more careful study and modelling 
which is beyond the scope of this paper. However, it is 
clear from some simple estimates, that this massive spiral 
galaxy is close enough to the HI complex to induce in it a 
gravitational instability via mechanisms such as those sug- 
gested by Icke (1985), Noguchi (1988) and Olson & Kwan 
(1990a, b). The mechanism proposed by Noguchi (1988) 
generates a central bar, which in turn causes gas to sink 



towards the center of the disk and induce its collapse. No 
observational evidence is seen, neither for a central bar- 
like structure in the HI cloud of SBS 0335-052, nor for a 
central density concentration. The mechanism by Olson 
& Kwan (1990a,b) works via a large tidal increase of the 
inelastic collision rate of individual gas clouds, and their 
merging and collapse. It is basically a stochastic process. 

The most probable mechanism is that proposed by Icke 
(1985) involving tidal acceleration to supersonic velocities 
of gas layers in the external parts of an HI disk, with sub- 
sequent dissipation of the kinetic energy and loss of dy- 
namical stability of the gas disk. It takes only a moderate 
tidal force, and acts from larger distances in comparison 
with the two other mechanisms. Adopting a radius R = 
16 kpc for the HI disk and a circular velocity V of 40 
kms~^ at R, a dynamical mass within i? of 6 x 10^ Mq 
for the SBS 0335-052 system, a total dynamical mass of 
10^2 Mq for NGC 1376, and a sound speed of the HI gas 
of 10 kms~^, we find that for pericenter distances smaller 
than ^ 412 kpc, the gas in SBS 0335 052 can be acceler- 
ated to supersonic velocities. This distance is significantly 
larger than the present projected distance of 150 kpc so 
that sufficiently strong tidal forces of NGC 1376 on the HI 
complex are expected. Such a tidal triggering mechanism 
may explain the fact that the two star-forming c;enters be- 
gan nearly simultaneously, and are aligned along the line 
joining the centers of the two objects, symmetrically on 
either side of the center of the SBS 0335 052 system, just 
as the tidal forces of the moon on the earth causes high 
ocean tides on diametrically opposite locations on earth. 

Alternatively, we may consider that we do not have a 
single huge HI complex, but two distinct smaller HI com- 
ponents interacting tidally with each other and triggering 
star formation at the location of the two HI peaks. In this 
case, the tidal force Fi exerted by each cloud is propor- 
tional to M (cloud) / df where di is the projected separa- 
tion between the two BCDs. On the other hand, the tidal 
force F2 exerted by NGC 1376 is proportional to M (NGC 
1376) / where ^2 is the projected separation between 
the galaxy and the SBS 0335-052 system. Adopting M 
(cloud) = 2 X 10^ Mq, M (NGC 1376) = lO^^ Mq, di = 
22 kpc and d2 — 150 kpc, we find Fi / F2 = 0.6, i.e the 
tidal forces in the two cases are comparable and we cannot 
say for sure which of the two scenarios discussed above pre- 
vails. In the case of two distinct HI components, we find 
that the Icke (1985) mechanism will apply for pericenter 
distances less than ~ 27 kpc. 

Clearly, detailed hydrodynamical calculations of the 
tidal interaction between two HI clouds and between a 
massive galaxy like NGC 1376 and a primordial gas cloud 
similar to that associated with SBS 0335 052, arc needed 
to understand better the triggering of star formation in 
these systems. 

5.3. Star formation threshold 

Star formation in disk galaxies is thought to occur above 
a critical threshold of the gas surface density above which 
clouds will become gravitationally unstable against col- 
lapse (Toomre 1964; Kcnnicutt 1989). This critical den- 
sity is proportional to the velocity dispersion of the gas 
disk and the epicyclic frequency, which is twice the cir- 
cular frequency for regions with solid body rotation. The 
validity of the Kennicutt (1989) HI column density thresh- 
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old for onset of star formation has been checked for BCDs 
(e.g., Taylor et al. 1994) as well as for LSB galaxies (van 
der Hulst et al. 1993). Alternatively, as first mentioned 
by Skillman (1987), star formation in dwarf irregulars oc- 
curs when the column densities reach values of order 5 x 
10^°cm~^, with this threshold increasing with decreasing 
metallicity (Franco & Cox 1986). 

The observed peak column densities in the SBS 0335- 
052 system coinciding with the optical counterparts are 
near the upper end of the range of those reported for 
HII galaxies by Taylor et al. (1993, 1995) — (0.8 - 11) x 
10^°cm~^ — when compared to data on objects mapped 
with comparable linear resolution, of order 5kpc. These 
values are also in the same range observed for LSB galax- 
ies, — (1.0 - 10.5) X 10^°cm~^ — objects in which no 
significant star formation takes place, and observed by van 
der Hulst et al. (1993) with a similar linear resolution. 

Given our modest spatial resolution, all we can do is 
make an order of magnitude estimate of the expected 
threshold value and compare that to the observed (beam- 
smeared) ones. To obtain an estimate we assume that both 
peaks fall within the solid body part of the rotation curve, 
at a rotational velocity of about 20 kms^^ and a radial 
distance of 6 kpc. With a typical gas velocity dispersion 
Gy of 6 kms~^, we obtain for the critical surface density 
at the radius of the current SF burst a value equal to only 
1.3 Mq pc-2 (equation 6 of Kennicutt 1989 with a = 2/3), 
in comparison with measured values of 8.0 Mq pc~^ and 
6.0 Mq pc~^. With a viewing angle of the HI disk of less 
than 77°, the gravitational instability condition would be 
satisfied. This appears to be the case since in the simple 
model of a circular flattened HI disk, the observed axial 
ratio would correspond to a viewing angle of 70°. 



5.4. The role of a dark matter halo 

It has been shown (e.g. de Blok et al. 1996) that LSB 
galaxies are dominated by dark matter over the entire area 
where HI and stellar emission is detected. This dark mat- 
ter dominance is one of the main factors which allows an 
LSB galaxy to be stable against various perturbations (see 
e.g., the study of the relative fragility of LSB and high- 
surface-brightness disks by Milios ct al. 1997). If in the 
SBS 0335-052 system we indeed observe the formation of 
a first generation of stars in a low density primordial gas 
cloud, then its dark matter halo should have played an es- 
pecially important role for stabilizing the gas against var- 
ious local and global instabilities and preventing it from 
gravitationally collapsing and forming stars. Presumably, 
a strong tidal perturbation due either to a close passage 
of a massive galaxy or to the interaction with another HI 
cloud can overcome the stabilizing infiuence of the dark 
matter halo and trigger a star burst. 

6. CONCLUSIONS 

From the observational data presented here, and the 
properties derived for the HI gas and based on the discus- 
sion above we arrive at the following conclusions: 

1. The optical dwarf galaxy SBS 0335-052 is associated 
with a huge HI cloud with overall size of 66 x 22 kpc, 
elongated in the EW direction. The cloud has a total 
HI mass of 1.68 x 10^ Mq. 

2. This HI cloud contains two prominent HI peaks, lo- 
cated nearly symmetrically to the East and West, 
and separated by about 22 kpc. They measure 7 x 
4 kpc and contain 0.79 xlO^ Mq and 0.89 xlO^ Mq, 
respectively, of HI. 
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3. The two HI density peaks are identified with two 
optical emission- line dwarf galaxies, the eastern one 
with SBS 0335-052 itself, and the western one with 
a new, one order of magnitude less luminous dwarf 
named SBS 0335-052W. This latter galaxy shows 
very similar properties, including an extremely low 
metallicity of the HII region gas, consistent with the 
amount produced during a first SF burst. Radial 
velocities of these galaxies, as measured from the 
emission lines of their HII regions are very close to 
those of the HI gas at the positions of their respective 
peaks. 

4. The velocity field of the HI cloud is rather disturbed, 
with presumably tidal tails protruding from both 
edges of the cloud. This suggests that what we are 
seeing is an interaction, possibly with the nearby 
galaxy NGC 1376. This large spiral galaxy, the dom- 
inant member of the galaxy group LGG 103, with a 
projected radial velocity of 120kms~^ larger than 
that of SBS 0335-052, is located at a projected dis- 
tance of 150 kpc. Consideration of the relevant pa- 
rameters of the SBS 0335-052 system and NGC 1376 
does not contradict the hypothesis that the forma- 
tion of the two young dwarfs from the huge low- 
surface density HI cloud was triggered by a tidal in- 
teraction with NGC 1376. 

5. Alternatively, we are dealing with two almost equal 
mass HI clouds, each measuring ~ 7 x 4 kpc, 
each containing an optical counterpart, and showing 
solid-body rotation with an amplitude of ~ 35-40 
kms~^. In this case, the triggering of the star for- 
mation in the two HI clouds is probably due to their 
mutual tidal interaction. 



6. Each component seems to require the presence of 
dark matter, and the system as a whole is domi- 
nated by a dark matter halo. The total dynamical 
mass of the system is estimated to be between 3.8 
X 10^ M© and 9 x 10^ Mq, or - 1.9 to 4.5 times 
larger than the total mass of gas and stars. 

7. If the viewing angle of the HI disk is less than 77°, 
the gas surface densities for both HI peaks are larger 

than the critical value determined from the param- 
eters of the rotation curve at those radii and are 
consistent with the onset of gravitational instability 
and subsequent star formation. 
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